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We analyze continuous-time quantum walks on necklace graphs - cyclical graphs consisting of 
many copies of a smaller graph (pearl). Using a Bloch-type ansatz for the eigenfunctions, we block- 
diagonalize the Hamiltonian, reducing the effective size of the problem to the size of a single pearl. 
We then present a general approach for showing that the mixing time scales (with growing size of 
the necklace) similarly to that of a simple walk on a cycle. Finally, we present results for mixing on 
several necklace graphs. 



I. INTRODUCTION 



o 

(N 

Q_j Classical random walks on form the basis for many successful physics-inspired algorithms. The evolution of prob- 

■ ability distributions according to simple update rules for probability spreading allows us to sample from thermal 
\ distributions (via the Metropolis algorithm [H, 0) or to look for ground-states of physical systems(with simulated 

■ annealing 0, |j] ) ■ The effectiveness of random- walk based algorithms can be characterized by its mixing time (how 
fast it approaches the stationary distribution), or by a hitting time (how fast it reaches a particular vertex). For 
example, the fast mixing of a random walk algorithm for sampling from the thermal distribution of the Ising model 
forms the basis of a fully polynomial randomized approximation scheme for the permanent of a matrix U 

Oh Thinking about how to utilize the probabilistic nature of quantum mechanics, instead of analyzing the diffusion 
of probabilities, we can ask what happens if we let the amplitudes in a system whose interactions respect some 
€ ' H^£j^ structure evolve according to the Schrodinger equation. The result of this way of thought are quantum walks 
[71, |8| , a useful tool in quantum computation. They bring new dynamics (different wavepacket spreading [ Qj) and 
qh| algorithmic applications (e.g. in searching for graph properties [l0| . graph traversal game evaluation as 
well as theoretical results (universality for computation |13|). We can define quantum walks in discrete time with an 
additional coin register, or in continuous time, with Hamiltonians which are adjacency matrices of graphs. In this 
. paper, we choose the latter approach. 

' The mixing of quantum walks has been previously inve stig ated for several types of graphs - e.g. on a chain Q, a 
\ 2D lattice [l3| , hypercubes [l^ and circulant graphs [IB, [i3l • this paper we focus on continuous quantum walks 
■r^lj- ■ on necklaces - cyclic graphs composed from many {K) copies of a subgraph of size M (pearls), as depicted in FigH] 
\ Our goal is to provide a simplified approach for finding their eigenvectors and eigenvalues, as well as for analyzing 
the mixing times for such walks. 

The motivation for analyzing this type of graph comes from Hamiltonian complexity [Tsj . Quantum computation in 
the usual circuit model |19| can be translated into a quantum walk in two ways. First, following Childs [l3|j evolving 
a wavepacket on a graph with many wires (representing basis states), connected according to the desired quantum 
circuit. Second, we can use Feynman's idea |20[ to view a computation as a "pointer" particle doing a quantum walk 
(hopping) in a "clock" register, while the computation gets done in a "data" register j2l| - |24| or particles holding the 
/\ \ working data hopping along a graph H^H^. In both cases, we need to look at transmission/reflection properties of 
the graphs, and their long-term dynamics. Specifically, we would like to know (and ensure) that a computation is 
done when we want it to be, not having the wavepacket localized (or spread) in undesired parts of the graph. This is 
why we focus on the mixing properties of quantum walks that arc underlying quantum computational models based 
on quantum walks, looking at their spectra in detail. Note that proofs of computational complexity for QMA-hard 
problems (e.g. HH, [13) ^-l^o involves investigating the (low- lying) spectrum of a quantum walk. The simplest graph 
involved in the Feynman-like models is a line or a cycle, and the dynamics for this quantum walk are well understood 
[3, [l^. We look at continuous-time quantum walks on necklace graphs, which appear in the analysis of quantum 
computational models [13, 113, Ull that generalize the Feynman approach. Necklace graphs could also be viewed as 
implementing dynamics for quantum walks on imperfect cycles. 

Utilizing the cyclic structure of the necklaces, we propose a Bloch-type ansatz for the eigenfunctions, allowing us 
to obtain several results. First, in Section [TT| we reduce the problem of finding the eigenvectors and eigenvalues of 
the quantum walk on a necklace of K pearls of size M to diagonalizing a M x M matrix K times (compared to full 
KM X KM diagonalization. Second, in Section [III[ we analyze average-time mixing for quantum walks on necklaces 
and find a general method for showing convergence to the limiting distribution. Finally, in Section IIVI we work 
out examples of quantum walks on particular necklaces, giving analytic (and numerical) results for the eigenvectors, 
eigenvalues and the scaling of the mixing time, concluding with open questions in Section IVl 
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FIG. 1: Examples of necklace graphs, with the pearls denoted by shaded regions. 
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FIG. 2: Using the plane wave-like ansatz Q in which neighboring pearls get a constant phase factor difference, the necklace 
Hamiltonian can be block-diagonalized on the pearls, acting independently on each pearl with the addition of a single link 
(carrying a phase factor) between its roots Xi and xm- 



II. FINDING EIGENVECTORS AND EIGENVALUES 



Consider a quantum system with a Hamiltonian H given by the adjacency matrix of a necklace-like structure. The 
simplest necklace is a cycle with K vertices. A general necklace is a collection of K pearls (small identical graphs 
with M nodes), connected into a cycle as in Fig. [TJ 

We label points in the j-th pearl xH'', with 1 < m < M. The endpoints of the j-th pearl (connected to the previous 
and following pearls) are x'f'^ and Let P be the adjacency matrix of a pearl. The Hamiltonian for the whole 
necklace is a sum of intra-pcarl terms and the connections between them: 

K K-l 



" " =1 



(I 



Our goal is to find the eigenvalues and eigenvectors of H. Because of the underlying cyclic structure of a general 
necklace graph with K pearls, we can assume that its eigenvectors will have a structure related to a plane wave on 
a cycle with K nodes. Let us then look at the if -node cycle first. There the Hamiltonian ([T]) has no P^'^^'s in it, 
allowing us to find the (plane- wave) eigenvectors of H°: 



K 

iPkj 



e 



X,) , (2) 



corresponding to eigenvalues parametrized by momenta pk ■ 

Xl=2cospl, Pl = ^^, (3) 

forfc==0,...,/\ - I. 

Consider now a general necklace with K pearls. Wc expect the eigenvectors of the necklace Hamiltonian ([T]) to have 
a form resembling also depending on the momenta pk dS]). Let us thus look for the eigenvectors of H in the form 



Hk) = ^f2e''^'\y[^^), (4) 
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where each 



M 



m— 1 

is a normalized vector with support only on the j-th pearl (the vertices x^"^ , . . . , x ) • Using dl]) and (g]) , we obtain 

1 



where the last two terms correspond to the amphtudes on the endpoints of the j-th pearl coming from the endpoints 
of the neighboring pearls. Notice that because of our parametrization the Hamiltonian is now block-diagonalized, 
acting in the same way on each pearl (see Figure [2|). When \ipk) is an eigenvector of H, we also have 



H\ipk) = Afe 



1 ^' 

\ipk) = Xk—;=y^ 



Using ^ and ([7|), finding the eigenvalues of H thus reduces to diagonalizing the M x M matrix 
where P is the adjacency matrix of a pearl, and 



(7) 



(8) 



ik = 







(9) 



has only two non-zero elements in the corners if a pearl has two distinct roots xi and xm- There is a special ease 
when a pearl is connected to the rest of the necklace through a single root vertex xi. There, the matrix Qf; has a 
single nonzero element and reads 



(10) 



Diagonalizing ^ gives us A/-dimensional vectors \yk)- For each fc = 0, . . . , /v — 1, there will be M of these, and we 
will label them \yk.n) with n = 1, . . . A/. The corresponding eigenvalues Xk.n of Yfc are also the eigenvalues of the full 
Hamiltonian H. Therefore, to find all the KM eigenvalues Xk.n of the necklace Hamiltonian with K pearls, we need 
to diagonalize the M x M matrix Yk jS]) for each A: = 0, . . . , A' — 1. To get the eigenvectors of H from the eigenvectors 

of Yk, we plug the coefficients ?/m'"^ of the vectors we just found into ^ and ([5]). 

In conclusion, the ansatz ([4]) simplifies the general problem of diagonalizing the KM x KM matrix H to diagonalizing 
an M X M matrix K times. This is useful especially when M is small and K is large. Our focus in what follows will 
be on mixing of continuous quantum walks on many-pearled (large-X) necklaces. 



III. QUANTUM WALKS AND MIXING 



A. Mixing in a time-averaged sense 



Time evolution according to the Schrodinger equation with a Hamiltonian that is an adjacency matrix of a graph 
produces a continuous time quantum walk. Let the eigenvectors of the system be {tpk) and the corresponding eigen- 
values Xk- When starting from an initial state \ipa), the probability of finding the "walker" at vertex \x) at time t 
(measuring position x) is 



p^>(t) = \{x\^m' 



-iXi,t , 



k 



x\xJjk){lJjk\VQ) 



(11) 
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The evolution is unitary, so it does not mix towards a time-independent stationary distribution like a classical Markov 
process. On the other hand, we can think about mixing for a quantum walk in a time-averaged sense, investigating a 
time-averaged probability distribution. It holds information about the probability of finding the system at a particular 
vertex at time <, chosen uniformly at random between and T (a chosen limiting time): 



p^r\T) ^ l-f p\r\T). (12) 

^ Jt=0 

This time-averaged probability has a well-defined T — > oo limit, which gives us the limiting probability distribution^ 
expressible using the eigenvectors of as: 

TT^" = hm p^r^T) (13) 



oo 



(14) 



k I \ •'^-^ J 

= E ('^i\^)i^\'^k){ipk\'Po){'Po\il^i), (16) 

where the final sum goes over pairs of equal eigenvalues. Note that for some quantum walks this limiting distribution 
can be dependent on the initial state (e.g. when we start in some eigenstate), so we will keep the superscript cpQ 
around. 

To determine how fast the time-averaged probability converges towards the limiting distribution, we need to bound 
the total distribution distance \\p'^°{T) — tt'-^"\\. Using (fT6)) . integrating an exponential and realizing that the terms 
summed over pairs of equal eigenvalues subtract out, we arrive at 



\\p'^^'{T)~n^"\\^J2\P-'-<'\ = E 



E/g-i(Afc-A,)T _ I 
V -H^fe - Xi)T 



(17) 



where the sum now goes over pairs of eigenvalues that are not equal. We can put an upper bound on this expression 
by a technique similar to [? ]. First, wc use |e~*(-*''''^'^')^ — 1| < 2 and move the absolute value inside the sums, to 
obtain 

The Cauchy-Schwartz inequality K-^/jx)! |(x|-!/'/c)| < \ (^l(V^iN)l^ + l(3;|''/'fe)l^j allows us to perform the sum over a;, 
resulting in 

Ib--(T) - vr-li < y 2 1(^.1^0)1 K^ol^OI^ ^^g) 

After another use of the Cauchy-Schwartz inequality on the terms involving |(/3o), realizing the expression is symmetric 
under exchange of k and we finally obtain 

lb--(r)-.-||< Y: (20) 

which corresponds to Lemma 4.3 of [? ]. It involves a sum of the inverse of eigenvalue differences. These terms can be 
large, but as T grows, the 1/T factor can bring the total variation difference to zero. It is our task now to investigate 
how fast this happens. We seek Tmix (the mixing time), for which 

||p'^«(r) -7r^«|| < e (21) 

would hold for all T > Tmix (e), given any precision parameter e. 
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B. Quantum Walk on a Cycle: The Limiting Distribution 

For our first example, we now follow [29| and compute the limiting distribution for the case of a walk on a cycle. 
Later, we will show that the time-averaged probability converges to it for times T ~ 0{K log A'), using a more general 
mixing result proved in Section [ill CI 

The eigenvalues and eigenvectors for the continuous-time quantum walk on a cycle arc given by ([2]) and We 
obtain the limiting distribution from (|16p by summing over the few nonzero terms. The sum over the equal eigenvalues 
splits into a sum over k — I and k + 1 = K (degenerate eigenvalue pairs). When the initial state \ipo) is concentrated 
at a vertex z, in the case of even K, the limiting distribution for the quantum walk on a cycle is 



-w-) = ^r^-]^' (22) 



where f^.z = for all pairs {x, z), with an exception for the two points x* = z and \x* ~ z\ = ^ , where its value is 
fx*.z = 1- For a cycle with an odd length K, we get 



= (23) 



with z defined in the same way as for even K , equal to zero for all pairs (x, z) except for x* = z, where fz^z = 1- 
The slight differences from a uniform distribution arise because not all of the eigenvalues are doubly degenerate. 

Proving that the time- averaged distribution converges towards the limiting distribution for T > 0{K log K) takes 
more work. We want to show that the total distribution distance \\p''''°{T) — Tr'-P°\\ goes to zero as T > 0{K log K). 
When computing p^°{T), the terms with Afc = Aj produce the limiting distribution tt*^" and are thus subtracted out. 
However, the terms left over (which were killed by the T — >■ cx) limit when computing tt) need to be carefully accounted 
for and govern the convergence. In ((20)) . we have a bound on the total distribution distance by a sum over pairs of 
inequal eigenvalues. We will upper bound this sum in Section IlII C[ using a general approach of lower bounding the 
terms |Afe — A^j^^ in (PP)) . This result is then applicable to several other walks on necklaces. 

C. A general approach to proofs of mixing 

The rate of convergence of the time-averaged distribution towards the limiting distribution is governed by a sum 
of \Xk.n — Xj,m\^^ over non-equal eigenvalues as in (|20p . We will now show a method for upper bounding it that will 
work in several cases. 

First, let us choose two particular sectors of eigenvalues, fixing n and m. It is often possible to bound the eigenvalue 
differences for this sector as 



|Afe,„ - Aj>i| > Cn,m IcOSPj - COSpfc| = 2Cn,r. 



sm 



7r(fc - i) 
K 



sm 



K 



(24) 



for some constant c„_m, where where pk = are the momenta, and the indices j, k run from to — 1, observing 



\k — j\ 7^ {0, yI- We will rewrite ([^^ using the substitution a = j + k and b = k — j. 

— yy, — , (25) 



V 1 ^ 

WA.,„ - 2c„,™ |sinf llsin^l 



where |6| < — 1 while 6 7^ 0, and |6| < a < 2{K — 1) — |6| while a ^ {0, K}. First, because of symmetry, we observe 
that it is enough to sum only over < & < ^^'^ multiply the resulting sum by 4. Second, it can only increase our 
upper bound if we count all a > 0, instead of having to take care with counting starting at The symmetry of the 
term involving a then again allows us to sum only up to a = and multiply the result by 4. Therefore, we obtain 

^ ' <T7^\}^T-^\\y:T—^\- (26) 





Recalling now that sin& > — , wc can deal with each sum as 



^ 1 K 

> T rT<— > -<— n— , 27 
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FIG. 3: A {K, c?)-comb graph is a ring of length Kd, with an extra vertex connected at every d-th node. 



thus finally expressing the sum in ()20p (note that we worked only in a single n, m sector) as 



E 



< 



2 KHn^^ 



r|Afc,„-Aj>„| TK I6c„ 



8c„ 



T 2 ' 



(28) 



where the extra factor K^^ comes from the term |('i/;fe|</'o)| , when we expect the initial state to have roughly equal 
overlap with all momentum states. According to (|20p and working this out for all sectors n, m, this suffices to show 
an upper bound on the mixing time (in the time-averaged sense) for this type of quantum walk, which grows with 
the system size a Tmixi^) < 0(e^^/l In^ A'). 

We will now show that for particular examples of walks on necklaces, the eigenvalues obey (|24p . and so that we can 
use the above approach for proving their convergence. The first example that we can deal with using this method 
is the cycle itself, where (|24p is an equality. Therefore, we have just shown that the time-averaged distribution 
(when starting from a single point) converges to the limiting distribution for a cycle of length K with mixing time 
Tmixi^) < O* (e~^if) (up to logarithmic factors). 



IV. EXAMPLES: QUANTUM WALKS ON COMB-LIKE NECKLACES 

We now look at a specific type of necklaces, which appear in the quantum- walk based model of computation (27| . 
These "combs" are constructed from a ring of length Kd by attaching an extra node (tooth) to the ring at every d-th 
vertex as in Fig. [31 The "pearl" in this comb- like graph has size d-l- 1, and there are K of them, so the total number 
of vertices in this graph is iV = K(d + 1). We will now analyze the spectra and mixing properties on the {K, d)-comb 
necklaces, showing their similarity to (and differences from) a simple cycle. 



A. The {K, l)-comb necklace 

This is the simplest of the graphs, with a pearl that has only two nodes (the base and the tooth): 



1 

1 



Because P has a single root, the matrix Qk needed to construct Yk ^ reads 



Therefore, the matrix Yi. IHl) is 



Its eigenvalues are 



and the corresponding eigenvectors are 



Y 



e'P" + e-'P" 




1 



X[^l^^ = COSpfe ± CQS^Pk, 



,(fe.±)\ = 



1 + K,± 



Xk.± 

1 



(29) 



(30) 



(31) 



(32) 



(33) 
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where the upper vector component corresponds to the base (and the lower component to the tooth) of the comb. 
According to this gives us the eigenvectors of H as 



IV' 



ik,±)\ 



1 



K 



1 + 



1 



(34) 



with pk = 2^ for k = 0, . . . , K —1. The eigenvalues ((5^ of the Hamiltonian are symmetrically distributed around 
zero, and each of them is also doubled if K is even. Note that two of the eigenvectors are zero on every other base 
and tooth, and correspond to eigenvalues ±1. 

The limiting distribution is analyzed in Appendix [Xj We find that for large K, the limiting distribution (when 
starting from a base vertex) is on base vertices and on teeth, with corrections for the initial vertex and 

the vertex across from it. We now prove convergence to the time- averaged limiting distribution, showing that the 
(time-averaged) mixing on this densest comb is no different than the one we saw for a cycle. We will upper bound the 

sum in ([20)) by the method in Section UllCi dividing the eigenvalues into 4 regions, -t-+, H — , — +, , corresponding 

to choices of n and m. In the H — and — h regions, we have |Aj.o — Afe.ij > const., so the inverse of such terms does 

not govern the scaling in (|20p . The important region combinations must then be and , where a few lines of 

algebra give us 



|A 



Afc,±| 



ipj - cospfcl . \Xj^± + Afc.il ^ \/2 - 1 

V2 



+ COS^Pfc 



iPj ~ COSpfcl 



(35) 



as the eigenvalues are well bounded away from zero. Armed with this inequality, and the fact that the overlap of a 
single-starting- vertex initial state with the eigenvectors scales as we can now use the result of Section UlI CI This 

gives us an upper bound on the mixing time for the (K, l)-comb necklace, scaling with K as T'^^^^(e) < O* (e^^A'), 
i.e. the same as for a cycle with no teeth, up to logarithmic factors. 



The {K, 2)-comb 



The next example is the (K, 2) comb. It has a tooth (extra vertex) at every second node of the basic loop, so its 
pearl P has three vertices. We label the base of the tooth as vertex 1, and the top of the tooth as vertex 2, giving: 



P = 



1 1 

1 
1 



(36) 



Following the procedure of Section [TTl we need to find the eigenvalues and eigenvectors of the matrices y'^'^^ con- 
structed as in JHl): 



After some algebra, we find that its three eigenvalues are 



'o 


1 


1 " 










f,-ipk - 







1 


e 2 2 cos ^ 


1 








+ 













1 








1 




















6^2™^,?"= 









A 



fc.O 



0, 



A 



(A'2) 
fe,± 



±^3 -f- 2cospfe, 



with the corresponding eigenvectors 

1 



\y 



(fc,0)\ 



^(3 + 2cospfc) 



2e~ 



2 COS ^ 



-1 



,(fc,±)\ 



v/2(3 + 2cospfc) 



±y/3 -f 2cospfc 
1 



2e- 



Pk, 



(37) 



(38) 



(39) 



To construct the eigenvectors of the Hamiltonian H, we use ([5^ in equation (U). 

Let us now look for a lower bound on the gap between eigenvalues. When we choose two eigenvalues from different 
sectors (0, + or — ), the differences between them are always larger than 1. The only interesting cases are thus the 
++ and choices of eigenvalue pairs. There we find 



|Aj,±-Afc,±|= ^^3 + 2cospfc - ^3 + 2cospj 



2 (cospfc — cos pj] 



1 



^3 -I- 2cospfc + ^3 + 2 cos Pj 2 



> t: (cospfe - cospj) , 



(40) 
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FIG. 4: A log-log plot of the smallest nonzero eigenvalue differences for {K, d)-comb graphs (see Figure [3| and a cycle. A 
straight-line fit through the datapoins indicates a K^^ scaling with the growing number of pearls. 



This lower bound on the nonzero eigenvalue gaps allows us to use the results of Section IIII CI and prove the upper 
bound T^,^^\e) < O* (e'^K) on the mixing time for the {K, 2)-comb. This is onee again the same upper bound we 
found for the cycle and the {K, l)-comb in Section HV Al 



C. The (Jf, d)-combs 

In the last example we want to show that comb-like necklaces with d vertices between teeth (see Figure mix 
similarly to cycles. We dealt with the most non-cycle-like examples in the previous Sections, and now we will 
numerically look at combs with general spacing d. The results for the smallest nonzero eigenvalue differences for 
d < 15 combs are shown in Figure HI In a log-log plot of the smallest eigenvalue difference vs. the number of pearls 
K (for various values of d), we see the characteristic K^^ scaling. Thus, the numerics imply that the scaling of the 

mixing time gets no worse than T^^^\e) = O (e~^if ^ log^ iiT) . However, it is likely that the eigenvalue differences 
also obey the cos-like scaling as we have seen for d = 1, 2. If we could show this, we would again prove that the 

mixing time scales with K as Tj^^'^\e) — O* (e^^A'). 



CONCLUSIONS 



The goal of this paper was to utilize the cyclical repetitive structure of necklace-like graphs, providing a general 
method for analyzing the eigenvectors and eigenvalues of continuous-time quantum walks on such graphs. Using a 
Bloch-theorem-like ansatz, we block-diagonalized the Hamiltonian, decreasing the effective size of the problem from 
KM to M, where M is the size of a pearl and K is the number of pearls in the necklace. Next, we wanted to 
investigate the mixing times (for approaching a limiting distribution in a time-averaged sense) for these quantum 
walks. In Section [ill CI we have shown that proving a cos- like lower bound (j24p on (non-equal) eigenvalue differences 



results in a mixing time T^ixi^) ^ O* (e^^A) for a graph with K pearls, which is the same as for a cycle with K 
nodes. Note though, that the prefactor in the mixing time can depend on the size of the pearl M . Finally, in Section 
IIVI we exhibited the bound (and thus the cycle-like mixing time) for two necklace-like graphs. These graphs 
appear in the models of quantum computation [24 . \2l\ that extend the Feynman-computer with the so-called railroad 
switches, and finding the polynomial-time (in K) scaling of the mixing-time is required for showing their effectiveness. 
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Appendix A: The limiting distribution for a {K, l)-comb necklace 

We now analyze the limiting distribution for the continuous-time quantum walk on a (K, l)-comb in Fig.O starting 
from a single (base) vertex z. First, we do our work analytically, and end with a few high- A' numerical approximations. 
The conclusion is that the distribution is flat except for a few points very close to the two special positions x* = z 
and \x* — z| = -y. 
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First, wc look at the time-averaged limiting distribution for going from base Xb to base Zb (which thanks to the 
identity (|A3p turns out to be the same as for going from tooth Xt to tooth Zt), obtaining 



where the last term involving K/2 occurs only for even K. Using the identities 



and denoting 



(Al) 



Afe,+Afe,_ ^ -1 (A2) 
11 A'^ 



1 1 

K 2(l + cos2pfc)' ^^^^ 
X A.2(l + cos2p,)' ' ^""^^ 



f for odd A' , 

C=\ (A6) 
I 2lf ^"-"^ even isT, 

with fx,z defined in ([22| . we rewrite (|A1|) to finally obtain 

TT (Xfalzfa) = TT (Xt|zt) = -1 (1 - A - B(^l^) + C + /^,^) . (A7) 

Next, we compute the time- averaged limiting distribution for the "start from a base - go to a tooth" transition, 
using the identity 

if - y — = - — (A8) 

+ 2(l+COS^p.) 

and (IMl). We obtain 



, , , , I , / ^(^ + ^N-) + 4l?) for odd i^, 
[ j^{A + S(^i^) + for even K. 

Finally, let us look at the high K approximation (see also FiglS]). In A we replaced the sum by an integral and 
obtained A « At points z for which [x — z) £ {0, the expression B(x|2) equals A exactly, while it rapidly falls 
off to zero with growing distance of z from x (or ^ + x). The limiting distribution when starting from a base point 
for even number of pearls is thus well approximated by a flat distribution with ~ on bases and — 5^ 

on teeth, with the exception of the starting pearl and the pearl exactly opposite to it receiving — and 

"^^^ — respectively. It is very similar for odd K, except that we do not have the special case of the opposite 
pearl. For an example of the limiting distribution with odd and even K, see Fig. [5] 
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FIG. 5; The limiting distribution p3|l on bases and teeth of a (K, l)-comb when starting from the base of the 50— th pearl for 
odd K = 201 (left) and even K = 200 (right). 



